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Abstract 

The characterization of the notions of complex and detailed balancing for mass action kinetics chemical 
reaction networks is revisited from the perspective of algebraic graph theory, in particular Kirchhoff’s 
Matrix Tree theorem for directed weighted graphs. This yields an elucidation of previously obtained 
results, in particular with respect to the Wegscheider conditions, and a new necessary and sufficient 
condition for complex balancing, which can be verified constructively. 


1 Introduction 

The notion of complex balancing of mass action kinetics chemical reaction networks, generalizing the classical 
notion of detailed balancing, dates back at least to the origin of chemical reaction network (CRN) theory; 
see especially [12, Til IS]- The assumption of existence of a complex-balanced equilibrium has powerful 
consequences for the dynamical behavior, precluding multi-stability and oscillations. In this note we will 
revisit the notion of complex balancing by a systematic use of notions and results from algebraic graph 
theory, in particular the Laplacian matrix and Kirchhoff’s Matrix Tree theorem (see also mm for other 
uses of this theorem in chemical reaction dynamics). This will result in a constructive necessary and sufficient 
condition for complex balancing. Furthermore, motivated by recent work in [4] expanding on [7], we will 
provide a new perspective and results on the Wegscheider conditions for detailed balancing and formal 
balancing as introduced in [3]. 

The structure of this note is as follows. Section 2 gives a brief recap of the basic framework of CRN 
theory from an algebraic graph theory perspective, based on mmm- Section 3 introduces Kirchhoff’s 
Matrix Tree theorem and shows how the application of this theorem leads to an improved, and more directly 
verifiable, condition for complex balancing as compared to CM- Section 4 relates Kirchhoff’s Matrix Tree 
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theorem to the notion of formal balancing, cf. [2] and [7J, and shows how this leads to an insightful graph- 
theoretic proof of the result obtained in [2] that complex balancing together with formal balancing implies 
detailed balancing and conversely. 

Notation: The space of n-dimensional real vectors consisting of all strictly positive entries is denoted by 
R” . The mapping Ln : R” — > R n , x H)• Lni, is the elementwise logarithm, and is defined as the mapping 
whose z-th component is given by ln(xj). Similarly, Exp : R n —> R” is the mapping whose z-th component 
is given by exp Xi. Furthermore, for two vectors x, y € R™ we let | denote the vector in R" with z-th 
component . Finally, l n denotes the n-dimensional vector with all entries equal to 1. 

Some graph-theoretic notions [3J: A directed graph@ Q with c vertices and r edges is characterized by its 
c x r incidence matrix , denoted by D. Each column of D corresponds to an edge of the graph, and contains 
exactly one element 1 at the position of the head vertex of this edge and exactly one —1 at the position 
of its tail vertex; all other elements are zero. Clearly, tjD = 0. The graph is connected if any vertex can 
be reached from any other vertex by following a sequence of edges; direction not taking into account. It 
holds that rank D = c — £, where i is the number of connected components of the graph. In particular, Q 
is connected if and only if ker D 1 2 3 = spanl c . The graph is strongly connected if any vertex can be reached 
from any other vertex, following a sequence of directed edges. A subgraph of Q is a directed graph whose 
vertex and edge set are subsets of the vertex and edge set of Q. A graph is acyclic (or, does not contain 
cycles) if and only if ker D = 0. A spanning tree of a directed graph Q is a connected, acyclic subgraph of Q 
that contains all vertices of Q. 

2 Recall of complex-balanced chemical reaction networks 

In this section, in order to set the stage, we will briefly recall the well-established framework of (isothermal) 
chemical reaction network (CRN) theory, originating in the work of Horn, Jackson and Feinberg in the 
1970s [12, [6] IU [9]. Consider a chemical reaction network with m chemical species (metabolites) with 
concentrations x £ R+, among which r chemical reactions take place. The left-hand sides of the chemical 
reactions are called substrate complexes and the right-hand sides the product complexes. To each chemical 
complex (substrate and/or product) of the reaction network one can associate a vertex of a graph, and to 
every reaction (from substrate to product complex) a directed edge (with tail vertex the substrate and head 
vertex the product complex). Let c be the total number of complexes involved in the chemical reaction 
network, then the resulting directed graph Q with c vertices and r edges is called the graph of complexes ^ , 
and is defined by its cx r incidence matrix D. Furthermore we define the mxc complex composition matriqj 
Z with non-negative integer elements expressing the composition of the complexes in terms of the chemical 
species: its k -th column expresses the composition of the fc-th complex. The dynamics of the chemical 

1 Sometimes called a multigraph since we allow for multiple edges between vertices. 

2 In the literature sometimes also referred to as reaction graphs. 

3 In [111 1121 . it is called complex matrix and is denoted by Y. 
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reaction network takes the well-known form 


x = Sv(x) = ZDv(x ), 


(1) 


where v(x) is the vector of reaction rates, and S = ZD the stoichiometric matrix. 

The most basic way to define v{x) is mass action kinetics. For example, the mass action kinetics reaction 
rate of the reaction X\ + 2 X 2 —>■ X 3 is given as v(x) = kx\x^ with k > 0 a reaction constant. In general, for 


a single reaction with substrate complex S specified by its corresponding column Z$ = 
the complex composition matrix Z , the mass action kinetics reaction rate is given by 


l T 


Zs\ 


Zsr. 


of 


1 rv.^S2 ' ' ' rp^STTl 

n,Uj 1 *^2 5 


which can be rewritten as fcexp(ZjLnx),® £ M™. Hence the reaction rates of the total reaction network 
are given by 


Vj(x) = kj exp(Zj.Lnx), j = 1, • • • ,r, 

where Sj is the substrate complex of the j-th reaction with reaction constant kj > 0. This yields the 
following compact description of the rate vector v(x). Define the r x c matrix K as the matrix whose (j, a)- 
th element equals kj if the er-th complex is the substrate complex for the j-th reaction, and zero otherwise. 
Then v{x) = ItExp (Z^Lnx), x £ MT, and the dynamics of the mass action reaction network takes the 
form 


x = ZDKExp (Z T Ln x), xe (2) 

It can be easily verified that the c x c matrix L := —DK has nonnegative diagonal elements and nonpositive 
off-diagonal elements. Moreover, since tjD = 0 also 1 = 0, i.e., the column sums of L are all zero. Hence 

L defines a weighted Laplacian mafriaQ for the graph of complexes Q. 

The aim of CRN theory, starting with mmm, is to analyze the dynamical properties of d^J), and 
in particular to derive conditions which ensure a dynamical behavior which is independent of the precise 
values of the reaction constants (which are often poorly known or varying). This has culminated in the 
deficiency-zero and deficiency-one theorems (see e.g. 0), while a somewhat complementary approach is 
based on the assumption of existence of a complex-balanced equilibrium 01m generalizing the classical 
notion of a detailed-balanced equilibrium. 

Definition 2.1. A chemical reaction network {I|) is called complex-balanceif there exists an equilibrium 
x* £ M™, called a complex-balanced equilibrium, satisfying 

Dv{x*) = -LExp (Z T Lnx*) = 0 (3) 

4 In [T2t 111] (minus) this matrix is called the kinetic matrix, and is denoted by A. 

s Note that in older references, e.g. mm, a reaction network (or mechanism) is called complex-balanced if there exists a 
complex-balanced equilibrium for all positive reaction constants. 
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Chemically (JSJ) means that at the complex-balanced equilibrium x* not only all the chemical species, 
but also the complexes remain constant; i.e., for each complex the total inflow (from the other complexes) 
equals the total outflow (to the other complexes). 

The assumption of complex balancing has been shown to have strong implications for the dynamical 
properties of fl2J); see in particular the classical papers mm®. As detailed in [18], expanding on EQ, 
these properties can be easily proved by defining the diagonal matrix 

S(x*) := diag (exp(ZfLnx*)). =lj jC , (4) 

and rewriting the dynamics © into the form 

x = — Z£(x*)Exp (Z T Ln (—)), C(x*) := LE(x*). (5) 

x* 

The key point is that since £(x*)Exp (Z T Ln (jy)) = 0 for x = x*, and Exp (Z r Ln (fr)) = l c , the trans¬ 
formed matrix £{x*) satisfies 

£(x*)l c = 0, l^£(s*) = 0, (6) 

and thus is a balanced Laplacian matrix (column sums and row sums are zero). Together with convexity of 
the exponential function this implies the following key fact. 

Proposition 2.2. 7 T £(x*)Fxp ( 7 ) > 0 for any 7 G R c , with equality if and only if D 1 7 = 0. 

First property which directly follows [T 8 | from Proposition 12.21 is the classical result [121 iTl . b] that all 
positive equilibria are in fact complex-balanced, and that given one complex-balanced equilibrium x* the 
set of all positive equilibria is given by 

£ := {x** G R+ | S T Lnx** = S T Lnx*} (7) 

Furthermore, using an elegant result from [ 8 ], there exists for every initial condition xq G RT a unique 
x** € £ such that x** — xq G im S. By using the Lyapunov function 

G(x) = x T Ln (^ 7 ) + 0** - x) T l m ( 8 ) 

Proposition 12.21 then implies that the vector of concentrations x(t) starting from xq will converge to x**; at 
least if the reaction network is persistent The chemical interpretation is that G is (up to a constant) the 
Gibbs’ free energy with gradient vector ^f(x) = Ln (^r) being the chemical potentials. See e.g. i f-T t [ 2 T, 22 j 
for further informatior]j. 

6 The reaction network is called persistent if for every xo £ R+ the w-limit set of the dynamics © does not intersect the 
boundary of R+. It is generally believed that most reaction networks are persistent. However, up to now this persistence 
conjecture has been only proved in special cases (cf. m, eg, Ei and the references quoted in there). 

7 The form © also provides a useful starting point for structure-preserving model reduction 12111181 1X7] . 
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3 A graph-theoretic characterization of complex-balancing 

In this section we will expand on earlier investigations to characterize the existence of a complex-balanced 
equilibrium (see in particular eh. a, and the references quoted therein), and derive a new necessary 
and sufficient condition for complex balancing which can be constructively verified. 

First note that by the definition of Ln : R+ —> M m the existence of an x* G such that LExp (Z T Lnx*) = 
0 (i.e., x* is a complex-balanced equilibrium) is equivalent to the existence of a vector p* € M m such that 

LExp {Z T p*) = 0, (9) 

or equivalently, Exp (Z T p*) G kerT. Furthermore, note that Exp (Z T p*) G Mf_. 

In case the graph Q is connected the kernel of L is 1-dimensional, and a vector p G M5_ (the closure 
of the positive orthant) with p G ker L can be computed by what is sometimes called Kirchhoff’s Matrix 
Tree f/ieorenj®[ which for our purposes can be summarized as follows. Denote the (i,j)-th cofactor of L by 
Cij = (—1 )* + - ? Mj j j, where Mi j is the determinant of the (i,j)-th minor of L, which is the matrix obtained 
from L by deleting its i-tli row and j-th column. Define the adjoint matrix adj(L) as the matrix with (i, j)-th 
element given by Cj t . It is well-known that 

L ■ adj(L) = (det L)I C = 0 (10) 

Furthermore, since l^L = 0 the sum of the rows of L is zero, and hence by the properties of the determinant 
function it directly follows that Cij does not depend on i; implying that Cij = pj, j = 1, ■ • • , c. Therefore 
by defining p := (pi, • • • , p c ) T , it follows from (flOl) that Lp = 0. Furthermore, cf. [3J Theorem 14 on p.58], 
Pi is equal to the sum of the products of weights of all the spanning trees of Q directed towards vertex i. 
In particular, it follows that pj > 0, j = 1, • • ■ , c. In fact, p 0 if and only if Q has a spanning tree. 
Furthermore, since for every vertex i there exists at least one spanning tree directed towards i if and only if 
the graph is strongly connected, we may conclude that p G if and only if the graph is strongly connected. 

Example 3.1. Consider the cyclic reaction network 


C 3 



8 This theorem goes back to the classical work of Kirchhoff on resistive electrical circuits m ; see 0 for a succinct treatment. 
Nice accounts of the Matrix Tree theorem in the context of chemical reaction networks can be found in [141110] , In [4] Kirchhoff’s 
Matrix Tree theorem is mentioned and exploited in the closely related, but different, context of investigating how far complex 
balancing is from detailed balancing; see the next section. 
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in three (unspecified) complexes C\,C2,C^. The Laplacian matrix is given as 


1c~^~ I _ Jc _1^"^" 

i \ j A/3 A/A/3 

L = —kf kf + k~2 — kf 

_ 1c _ 1c~^~ 1c~^~ I Jc 

Ay 3 Ay 2 Ay 3 | Ay 2 

By Kirchhoff’s Matrix Tree theorem the corresponding vector p satisfying Lp = 0 is given as 


P = 


P P ~b I P P ~P I P p 

Ay 2 A/3 A/ A/3 Ay ^ A/2 

Jc~^~ Jc~^~ I 1c~^~ If I Jc Jf 

Ay 1 Ay 3 | A/ Ay 2 Ay 2 Ay 3 5 

I 1c~^~ I 1c ^ 1c I 1c 1c 

Ay 1 Ay 2 | Ay 2 Ay 3 “| A/-J^ Ay 3 


where each term corresponds to one of the three weighted spanning trees pointed towards the three vertices. 


In case the graph Q is not connected the same analysis can be performed on any of its connected 
components. 

Remark 3.2. The existence (not the explicit construction) of p £ satisfying Lp = 0 already follows 
from the Perron-Frobenius theorem m, m Lemma V.2J; exploiting the fact that the off-diagonal elements 
of —L := DK are all nonnegatirJ®[ 

Returning to the existence of p* £ R m satisfying LExp (Z T p*) = 0 this implies the following. Let Qj. j = 
1, ■ • • ,£, be the connected components of the graph of complexes Q. For each connected component, define 
the vectors p 1 , ■ ■ ■ , p^ as above by Kirchhoff’s Matrix Tree theorem (i.e., as cofactors of the corresponding 
diagonal sub-blocks of L or as sums of products of weights along spanning trees). Define the total vector 
p as the stacked column vector p := co^p 1 ,--- ,//). Partition correspondingly the composition matrix Z 
as Z = [Z\ ■ ■ ■ Z(\. Then there exists p* £ M m satisfying LExp (Z T p*) = 0 if and only if each connected 
component Qj, j = 1, ■ ■ • , L, is strongly connected and 

Exp (Zjp*)=p jp i, f3j > 0. (11) 

This in its turn is equivalent to strong connectedness of each connected component Qj and the existence 
of constants (3) such that Zjp* = Ln pi + /?' 1, j = 1,-- - ,£. Furthermore, this is equivalent to strong 
connectedness of each connected component of Q, and 

Ln p £ im Z T + ker D T (12) 


Finally, (fT2l) is equivalent to 

D t Ln p £ im D J Z J = im S T 


(13) 


Summarizing we have obtained 

9 This implies that there exists a real number a such that — L + al m is a matrix with all elements nonnegative. Since the 
set of eigenvectors of — L and — L + aim are the same, and moreover by 1 T L = 0 there cannot exist a positive eigenvector of 
— L corresponding to a non-zero eigenvalue, the application of Perron-Frobenius to — L + aim yields the result; see |20l Lemma 
V.2] for details. 
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Theorem 3.3. The reaction network dynamics x = — .ZXExp (Z T Lnx) on the graph of complexes Q is 
complex-balanced if and only if each connected component of Q is strongly connected (or, equivalently, p £ 
M+j and [T3\) is satisfied, where the elements of the sub-vectors p> of p are obtained by Kirchhoff’s Matrix 
Tree theorem applied to L for each j-th connected component of Q. Furthermore, for a complex-balanced 
reaction network a balanced Laplacian matriC(x*) defined in |]5|) is given as 

C{x*) = Ldiag(pi,--- ,p c ) (14) 

Remark 3.4. The above theorem is a restatement of Theorem 3C in m; the main difference being that in 
m the positive vector p £ kerL remains unspecified, while in our case it is explicitly given by Kirchhoff’s 
Matrix Tree theorem. 

We directly obtain the following corollary stated before in m Eq. (3.21)]: 

Corollary 3.5. The reaction network dynamics x = — ZLExp (Z T Lnx) is complex-balanced if and only if 
p € and 

PT ' P 2 2 •'' ' Pc c = E (15) 

for all vectors a = (<ti,<T 2 ,--- , (J C ) T £ kerZ flimH. In particular, if ker Z PI imD = 0 (zero-deficiency 
mmw then x = — ZLExp (Z T Lnx) is complex-balanced. 

Proof. Ln p £ im Z T + ker D 1 if and only if <r T Ln p = 0 for all a £ (irn Z T + ker D 1 )- L = ker Z fl irn D, or 
equivalently 


0 = (Ti In pi 4-1- a c ln p c = ln pf 1 -\ -b ln pf c = • • • pff) 


for all <7 £ ker Z n im D. 


Example 3.6. Consider Example \3.1\ for the special case k 1 = k 2 = k 3 =0 (irreversible reactions). Then 
the vector p reduces to 


P = 


i U+ 3U+ 

^2 ^3 

JU+ 1 ,+ 

3U+ Z-+ 
^1 ^2 


The reaction network with complex composition matrix Z is complex-balanced if and only kf > 0, i = 1,2,3, 
and 


D T Lnp £ im D T Z T , 


D = 


-1 

1 

0 


0 1 

-1 0 

1 -1 


°It can be easily seen that the balanced Laplacian matrices C{x*) for different equilibria ** just differ from each other by a 


positive multiplicative constant for each connected component of Q. 
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This last condition can be further written out as 


In 

In 

In 



€ im D t Z t 


As a mathematical example, take the 

complexes X\ +X2,X2,2X\ +X 2 ). In 
0, > 0, and (k±) 2 = k^k(f. 


complex composition matrix Z = 


(corresponding to the 


1 0 2 
1 1 1 

this case the network is complex-balanced if and only if kf > 0, klf > 


4 Relation with the Wegscheider conditions and detailed balancing 


In this section we will relate the conditions for complex balancing as obtained in the previous section to 
’Wegscheider-type conditions’. This will also relate complex balancing to the classical concept of ’detailed 
balancing’. 

Throughout this section we will consider reversible chemical reaction networks, in which case the edges 
of Q come in pairs: if there is a directed edge from vertex i to j then there also is a directed edge from j to 
i (and in the case of multiple edges from i to j there are as many edges from i to j as edges from j to i). 
This means that the connected components of Q are always strongly connected, or equivalently, that p as 
obtained from Kirchhoff’s Matrix Tree theorem is in M+. 

Define the undirected graph Q as having the same vertices as Q but half its number of edges, by replacing 
every pair of oppositely directed edges of Q by one undirected edge of Q. Denote the number of edges of Q 
by r = Tj'r. Endow subsequently Q with an arbitrary orientation (all results in the sequel will be independent 
of this orientation), and denote the resulting incidence matrix by D. Clearly, after possible reordering of 
the edges, D is related to the incidence matrix D of Q as 


D = 


D 



(16) 


To the j-th edge of Q there now correspond two reaction constants k~j~, k- (the forward and reverse reaction 

_ fcf 

constants with respect to the chosen orientation of Q). Then define the equilibrium constants K^ q := -A=,j = 

3 

1, • ■ • , f, and the vector K eq := (K( q , ■ ■ ■ , Kf i ) T . 

Recall O] that the reaction network is called detailed-balance (0 if and only if it satisfies 


Ln K eq € im S T , 


(17) 


where S := 


ZD is the stoichiometric matrix of the reversible network with graph Q. This is equivalent to 


o\ lnid( q + • • • + <7flnRT r 5 q = 0 

11 This means, see e.g. I15lll2lf2l] , that there exists an equilibrium for which every forward reaction is balanced by its reverse 
reaction. 









(18) 


for all er = (<n, • • • , <r f ) such that o T S T = 0. Writing out K) q = -L this is seen to be equivalent to 

3 

= (*rr ■■■(** r 

for all o such that So = 0, known as the (generalized) Wegscheider conditions. 

Recently in [3] the notion of formally balanced was introduced, based on Feinberg’s circuit conditions in 
[7], and weakening the above Wegscheider conditions. In our set-up this notion is defined as follows. 

Definition 4.1. The reversible reaction network Q with incidence matrix D, and vector of equilibrium 
constants K eq is called formally balanced if 

Ln K eq € irn D T 

Since S = ZD ’formally balanced’ is trivially implied by ’detailed-balanced’, while if im5 T = irn D 1 
(zero-deficiency) the reverse holds. Furthermore, any reaction network with acyclic Q (and thus ker D = 0) 
is automatically formally balanced. 

As above, the notion of ’formally balanced’ is seen to be equivalent to 
o i lniL( q + • • • + <7 r lnAA q = 0 

for all o = (cti, • • • , Of) T such that o 1 D 1 = 0, which in turn is equivalent to 

(ktr---(ktr = (kfr---(kfr ( 19 ) 

for all o such that Do = 0 (that is, for all cycles o). We will refer to (fl9l) as the weak Wegscheider conditions. 
Note that the weak Wegscheider conditions only depend on the structure of the graph Q (i.e., its cycles) 
and the equilibrium constants, and not on the complex composition matrix Z as in the case of the ’strong’ 
Wegscheider conditions (fl8l) . 

Theorem 4.2. Consider a reversible chemical reaction network given by the graph Q with incidence matrix 
D, and with p determined by L. The following statements are equivalent 

1 . Ldiag (pi, - ■ ■ , p c ) is symmetric 

2. Ln (A" eq ) = D T Lnp 

3. Ln(iL eq ) € im D T (formally balanced) 

Proof. (1) (2) 

Let us first prove the equivalence between (1) and (2). Consider the i-th and j- th vertex of Q, and suppose 
that the orientation has been taken such that the a-th edge between i and j is such that i is the tail vertex 
and j is the head vertex. Then the (i, j)-th element of Ldiag (pi, • • • , p c ) is given by k~pj , while the (J, z)-th 
element equals kffpi. Symmetry of Ldiag (pi, • • • ,p c ) thus amounts to 

k «Pj = k aPi 
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for all pairs of vertices i,j. On the other hand, the a-th element of the vector D T Lnp is given by 
In pj - In pi 

while the a-th element of Ln ( K eq ) is given by 

Z-+ 

ln — = In kt — ln k~ 
u- a a 

ha 

Equality of ln pj — ln pi and ln £;+ — ln k~ is thus equivalent to 
In pj -)- ln k~ =\npi + ln 

which in its turn is equivalent to k+pi = k~pj as above. 

(2) ^ (3) 

Obviously (2) implies (3). For the reverse implication, consider any pair of vertices linked by an edge of the 
graph Q. Depending on the orientation of Q refer to one vertex as the tail vertex t and the other vertex as 
the head vertex h. Refer to the positive reaction constant from t to h by k + and to the negative reaction 
constant by k~. Now consider a spanning tree directed towards t with product of weights denoted by Tt. In 
case the edge between t and h is part of this spanning tree then it follows that by reversing the orientation of 
this edge it defines a spanning tree directed towards to h with product of weights denoted by r^. It follows 
directly that 

pTt = T h (20) 

In case the edge between t and h is not part of this spanning tree then divide the edges of the spanning tree 
into two sets; the set E± containing the edges of the part of the spanning tree from t to h (containing say £ 
edges) and the set E 2 containing the remaining edges of the spanning tree. Observe that E\ together with 
the edge from t to h forms a cycle of the graph Q. Since Ln ( K eq ) e im D 1 it follows that for the reaction 
constants along this cycle (choosing an appropriate orientation), 

k + ■ k^ ■ ■ ■ kf = k~ ■ ki ■ ■ ■ kj. (21) 

Now within the spanning tree directed towards t, if the orientation of each of the edges of E\ is reversed, we 
obtain another spanning tree directed towards h with product of weights denoted again by Th ■ By using (12TT) 
it is readily verified that also in this case we obtain the same relation (12011 . Summing up over all spanning 
trees we thus obtain the equality 

pft = Ph, (22) 

which can be equivalently written as ln p- = lnp/j — lnp^. Doing this for all adjacent vertices t and h this 
exactly amounts to the required equality Ln(RT eq ) = D T Lnp. B 
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Example 4.3. Consider again the reaction network described in Examvle \3.1\ (without specifying the com¬ 
plexes Ci,C 2 ,Cs). The transformed Laplacian matrix is computed as 



A - ! - 1 A _A _ 

/ly I M '^ M ^ a 3 


k+k++k-k++h-k 2 0 0 

c = 

— k+ —k-2 


0 k+k^-\-k+k 2 +k 2 k% 0 


— ^3 — ^2 ^2 +^ 3 ” 


0 0 ifef k^ -\~k~2 k% +fc k% 


(k+ +fc+) (fc+ k++k~ k ++ k 

ik 2 ) — k 1 (k+ k^ -\-k+ k 2 ~\~k 2 k 3 ) — k^ (k^ k 2 -\-k 2 k^ -\-k 1 k 3 ) 


-kf(kfkf-\-k- L fc 2 ) ( k 1 +kf )(kf kf +kf k 2 +fc 2 ^3 ) —& 2 +k 1 fc 3 ) 

— k 3 (kf kf-\-k 3 kf+k^k 2 ) — (^1 ^3 +^1" k^ +k^k 3 ) (k 2 +kf)(kf kf+kf k 3 +fcj” k 3 ) 

which is symmetric if and only if 


h + h + h + — h~h~h~ 

Ay (V2 A/2 Ay Ay 2 Ay 2 


(23) 


On the other hand, Ln K eq G L> T amounts to 


r A*”*" 1 
In p- 

k i 


-1 

0 

1 

fc) - 





In t 2 - 

kr, 

A i 

ln pr 

L k 3 _ 

G im 

1 

0 

-1 

1 

0 

-1 


fc-T fc-r I" 

which reduces to ln -4- + In + ln -ir = 0, and hence to the same condition 

«1 fc 3 


. T/ims f/ie reaction network 


is formally balanced if and only if K23\l holds. 


Now let us relate all this to the necessary and sufficient conditions for complex balancing obtained before, 
cf. (fTHD . Note that by (fTUD 


D T Lnp G im D T Z T ^ D T Lnp G im D‘ Z 


T ryT 


Hence a reversible reaction network is complex-balanced if and only if D 1 Ln p G im D 1 Z 1 = imS 2 
We directly obtain the following corollary proved by other methods 


D r Ln p € 
in S 


Corollary 4.4. A reversible reaction network is detailed-balanced if and only if it is formally balanced as 
well as complex-balanced. 


Proof. We have seen before that ’detailed-balanced’ implies ’complex-balanced’ as well as ’formally bal¬ 
anced’. For the converse we note that formally balanced implies that Ln(iL eq ) G im D T . Hence by Theorem 
l4~2l Ln (K eq ) = D T Lnp. Since furthermore the network is complex-balanced D T Lnp G imS^. Hence 
Ln(iL cq ) = D r Lnp G imS T , i.e., the reaction network is detailed-balanced. ■ 

In case the reversible reaction network is formally balanced the symmetric matrix C{x*) = Ldiag {p \, • • • , p c ) 
can be written as 


Ldiag (pi, ■ ■ ■ ,p c ) = DKD 1 

12 The result is formulated in [4] as the equivalence between ’complex-balanced’ and ’detailed-balanced’ under the assumption 
of ’formally balanced’. 
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where 1C is the r x r diagonal matrix, with a-th diagonal element given by := k+Pj = k~pi where 
the a-th edge of Q corresponds to the reversible reaction between the z-th and the j-th complex. For the 
interpretation of the positive constants n a as conductances of the reversible reactions please refer to [5l [22]. 

If additionally the formally balanced reaction network is complex-balanced (and thus, cf. Corollary 14.41 
detailed-balanced), then its dynamics thus takes the form 

x = -ZD/CD T Exp(Z T Ln 

In this case, see e.g. caEu, all equilibria x** are in fact detailed-balanced equilibria, that is, D T Z T Lnx** = 
Ln (K eq ) (= D T Lnp). 

5 Conclusions 

By a systematic use of notions from algebraic graph theory, in particular the Laplacian matrix and Kirch- 
hoff’s Matrix Tree theorem, previously derived results on complex, detailed and formal balancing have been 
proved in a simple manner. Furthermore, it has resulted in a new necessary and sufficient condition for 
complex balancing, which can be verified constructively. 

The results obtained in this note can be immediately extended to mass action kinetic reaction networks 
with constant inflows and mass action outflow exploiting the classical idea of adding a ’zero complex’; see 
[9] and [23] for further details. 

Current research is concerned with the application of the developed framework to questions of occurrence 
of multi-stability and structure-preserving model reduction; see for the latter also [ZD EHJ El- 
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